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We derive an expression for the Casimir force between slabs with arbitrary dielectric properties 
characterized by their reflection coefficients. The formalism presented here is applicable to media 
with a local or a non-local dielectric response, an infinite or a finite width, inhomogeneous dissipative, 
etc. Our results reduce to the Lifshitz formula for the force between semi-infinite dielectric slabs by 
replacing the reflection coefficients by the Fresnel amplitudes. PACS: 12.20.D.s,03.70.+k,77.55.-|-f, 
Q ■ 78.67.-n 
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I. INTRODUCTION 



Even though Casimir predicted in 1948 an attractive force between perfectly conducting plates placed in quantum 
vacuum, [jlj it is only in recent years that experimental studies of Casimir forces have reached the necessary accuracy to 
test in detail theoretical predictions. The first measurements were done by Derjaguin etal. § in 1951 using dielectric 
materials. In the following decades, a number of experiments to measure Casimir interactions between dielectric or 
conducting materials were performed, however involving large relative errors in the measured forces |j. It was until 
1997 that Lamoreaux performed measurements with a precision of the order of 5 % Q by using an electromechanical 
system based on a torsion balance. Other experiments were made taking advantage of the sensitivity of atomic force 
microscopes achieving precisions close to 1% HJtJ. Additional measurements have been made by Chan et al. || using 
a micro torsional balance. This experiment is representative of the effects that Casimir forces have in micromechanical 
systems as was theoretically shown by Serry and Maclay j^] . Applications to nanostructures have also been considered 
[p0| |l2| . This has boosted investigations in which the detailed properties of the materials such as absorptivity, rugosity, 
or finite temperature effects are taken into account in the theoretical calculations of the Casimir forces fig] ] . 

The standard approach to study vacuum forces between imperfect conductors is the macroscopic theory proposed 
by Lifshitz in 1956 for semi-infinite dielectric materials. Jl4|] In this theory, the dissipative effects associated to the 
O^' radiation reaction of the elementary atomic dipoles composing the dielectric is balanced by the fluctuating vacuum 
field in accordance with the fluctuation-dissipation theorem. 

For a configuration of two semi- infinite slabs, with dielectric permitivities e\ and e 2 , separated by a gap of width L 
and permitivity 63 , the Lifshitz formula for force per unit area is 
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and 



Gl( ^) = C3*l+6ip6 3a2+ea p x e2ip ^ L/c (2) 
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Ga (£,p) = Sl+PS2+P x e^v^Llc. 1; (3) 
si -ps 2 -p 

where £ = iu is an imaginary frequency, p and Si t2 are defined in terms of the momenta parallel and perpendicular 
to the slabs k, and Kf = k 2 + £i(i£)£ 2 / c 2 , respectively: K 3 — y^^p/c and K 2 2 — e.z^ 2 s\ 2 /c 2 . Lifshitz theory 
has been succesfully employed in a number of experimental situations, and it yields the Casimir force for perfect 
conductors. However, as pointed out by Barash and Ginzburg, Ka | it is not clear how to generalize the theory to more 
complex problems, such as nonplanar surfaces, multilayer systems, anisotropic media, etc. Thus, they proposed an 
alternative approach to Lifshitz formula. As a dissipative system does not posses well defined natural frequencies of 
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oscillation, they introduced an auxiliary system in which the dielectric permitivity depended only parametrically on 
the frequency. This procedure enabled them to calculate the free energy of the field as a sum over allowed states of 
the system. In turn, the force per unit area was obtained as the derivative of the free energy. Kats p2| re-elaborated 
the formalism of Barash and Ginzburg by writing the dispersion law for surface electromagnetic waves in terms of 
the reflection amplitudes r s and r p of the media for s and p polarizations. Approximating the reflection coefficients 
in terms of the frequency and wavevector dependent surface impedance Z(u>, Q) he was able to obtain approximate 
non-local corrections to the Casimir force for good conductors. Kats remarked that dielectrics materials required 
alternative formulations, as their reflection coefficients cannot be expressed merely in terms of surface impedances. 
This is not necessarily correct, as an exact relation between surface impedance and reflection coefficients may indeed 
be introduced for arbitrary systems |25| . Jaekel and Reynaud ]2(| also rederived Lifshitz formula in terms of reflection 
coefficients for partially transmitting mirrors. Their expression reduces to that obtained by Barash and Ginzburg. 
p5| However, their derivation is not valid when dissipation is included. 



Other approaches have been used to study the vacuum fluctuations in the presence of dielectrics 16-|lq|. More 
recently, the problem of quantization in absorbing media and its applications to Casimir forces has been considered by 
Kupiszewska [ ^9|j2C| l and also by Matloob [^lj in the one-dimensional case. Interestingly enough, when temperature 
effects are neglected the expression for the Casimir force in absorbing and non-absorbing materials has the same 
functional form. This fact suggests that it is possible to obtain the Casimir force between two dispersive and absorbing 
slabs without the need of quantizing in an absorbing medium as was shown by Reynaud et al. ||(| using a scattering 
matrix formalism. 

Within the framework of the above discussion, it seems valuable to present an alternative, very simple derivation 
of the Casimir force, valid for materials with arbitrary dielectric properties. This is the purpose of the present paper. 

II. FORMALISM. 

Consider two slabs a = 1,2 parallel to the xy plane within free space and separated by a distance L along the 
z-direction, with inner boundaries at z\ = and Z2 = L. We assume that the slabs are non-chiral, translational 
invariant and isotropic within the xy plane, but otherwise they may be arbitrary; they could be identical to each 
other or different, they might have a local or a non-local dielectric response, an infinite or a finite width, they may 
be opaque or transparent, dissipative, inhomogeneous, etc. We want to describe the electromagnetic field only within 
vacuum, so we hide all the details of the field-matter interaction within the slabs in their reflection amplitudes r s a 
and r v a (a = 1,2). The reflection coefficients are determined by the generalized surface impedances Z" (a — s,p) 
through 



a Z2 + z«- 

Here Z£ is defined as the quotient E^/H^ of the components E"^ and of the a-polarized electric and magnetic 
fields evaluated at the a-th interface for outgoing boundary conditions beyond z a , taken along appropriately chosen 
directions parallel to the surface, Zq = q/k and Zq — k/q are the surface impedances of vacuum, q± = (Q,±fc) 
are the vacuum wavevectors with projection Q parallel to the surface and components ±fc normal to the surface, 
q = \Jq^ = u>/c, where uj is the frequency and c the speed of light. The sign of k is chosen so that <f+ propagates (or 

decays) as z increases. Upon each reflection, Q and u> are conserved while the sign of ±fe is reversed. Notice that Eq. 
(m is exact and that no approximation is involved by the use of our generalized surface impedances, unlike other works 
|p!3| , p4| that use an inappropiate definition of surface impedance. It should be noticed that for local homogeneous 
semiinfinite media, Z 8 a = q/k a and Z p — k a /(e a q), where k a is the component of the wavevector normal to the surface 
within medium a with local dielectric response e a (w), and Eq. (Q) yields the well known Fresnel amplitudes. However, 
Eq. ([I]) is much more general. [f25f 

The density of states within vacuum may be obtained from the Green's functions of the system. To this end, 
we study first the case of s-polarized waves choosing x — z as the plane of incidence. With that choice, E = 
(0,^,0), B = (B x ,0,B z ), q± = (Q,0,±k) and the boundary conditions for E become iqE y (Q + ) = -Zfd z E y (0 + ), 
and iqE y (L~) = Z|<9 2 _E a (X~), with d z = d/dz. The electric Green's function is 

«M=», (5) 



where z < and z > are the smaller and larger of z and z' respectively, 
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and 



E>(z) 



+ r s -ik(z-L) 
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are solutions of the scalar ID wave equation with wavenumber k obeying the appropriate boundary conditions at 
z = + and L~ respectively and W is their Wronskian. Analogously, B obeys iqB x (0 + ) = — (Z^) 2 d z B x {Q + )/Z{ 1 
and iqB x (L~) = (ZJ,) 2 d z B x (L~) /Zf, so that the magnetic Green's function is obtained by replacing E y — ► B x and 
r s a — ► — r* in Eqs. We do not consider B z separately, as it is simply proportional to E y . 

For each Q, the local density of states per unit k 2 is given by Effl 



so that by substituting Eqs. 



Ph (z) = -^Im (Gf 2 (z, z) + G% (z, z)) , (fc = + z0+) 
-(0) and its magnetic analogues we obtain 
1 + r{r s 2 e 2fkL 
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independent of z. The density corresponding to p polarization may be derived similarly, and is simply given by 
Eq. @ after replacing all the superscripts s — > p. Finally, the total density of states is p^ = p s k - 2 + p\i ■ 

A photon in a state characterized by a,Q and k 2 has momentum zLTik and moves with velocity ztck/q along the 
2 direction, so that its contribution to the momentum flux is tick 2 /q. Multiplying this by the photon occupation 
number, integrating over k 2 with the weight function p^ 2 and adding the contributions from all values of a — s,p and 
Q, with the usual replacement ]\2q . . . — > A/2it J QdQ . . ., we obtain the momentum flux from the vacuum gap into 
slab 2. There is a similar contribution coming from the semiinfinite vacuum on the other side of the slab, obtained 
by substituting r% — > above and reversing the flux direction z — > — z. The total force per unit area is obtained by 
subtracting the contributions from either side Eq] 
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The integral over k runs from iQ to and then to oo, so that q remains real and positive. It is easy to show that for 
perfect mirrors (r° — ±1) Eq. ( |Tc| ) yields the expected Casimir force. Furthermore, simple substitution of the Fresnel 
amplitudes 



s k k a 

r = 

k + k a 

and manipulation of the integration contours in Eq. ([!(]) leads to the formula of Lifshitz (Eq.(ffl)). 
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III. CONCLUSIONS 



We have derived a general expression for Casimir forces between slabs with arbitrary dielectric properties charac- 
terized by the reflection coefficients of the material. This procedure avoids complications related to the quantization 
of the electromagnetic field in dispersive media. The expression we obtain for the Casimir force is convenient for 
calculations since the reflection coefficients can be obtained straightforwardly in theoretical computations or through 
experimental studies. Our approach is based on the exact definitions of surface impedance and yields the Lifshitz 
formula for semi-infinite slabs when the reflection coefficients are replaced by the Fresnel amplitudes. This contradicts 
the results obtained by Mostepanenko and Trunov [^3) and also by Bezerra et al. p4j that claim that the use of surface 
impedances is only an approximation valid for small transverse wave vector Q. However, this conclusion arises from 
considering an approximate expression of the surface impedance. In our work this limitation is not present, since we 
do take into account that the surface impedances for the p and s polarized waves are different. Therefore, the results 
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that we obtain are valid for any value of the wave vector Q. The expresion for the Casimir force in Eq. (10) can 
be used to calculate accurately the force between non-homogeneous systems. In fact, current experimental setups for 
measuring Casimir forces consist of two metallic surfaces with a high reflectivity in the frequency range of interest 
lo ~ c/a, coated by thin (~ 8 nm) Au/Pd layer to avoid Al oxidation. The analysis of the experimental data is per- 
formed by giving arguments on the transparency of the Au/Pd layer, which allows the use of the Lifshitz expression 
for semi-infinite homogeneous media with semi-empirical corrections for the presence of that layer. However, within 
our formalism, such assumptions are unecessary. In a future work, we will present exact results for Casimir forces 
between heterostructured media, generalizing our previous one-dimensional analysis [ fTo| |. 

Acknowledgments: This work was partially supported by NASA Breakthrough Propulsion Physics Project, and by 
D G APA-UN AM Project IN-110999. 



[1] H. B. G. Casimir, Proc. Kon. Ned. Akad. Wet. 51, (1948) 793. 
[2] B. V. Derjaguin, I. I. Abrikosova, Vestnik Akad. Nauk. SSSR, 6 (1951) 125. 
[3] P. W. Milonni and Mei-Li Shin, Contemporary Physics, 33, (1992) 313. 
[4] S. K. Lamoreaux, Phys. Rev. Lett. 78, (1997) 5. 

[5] H. B. Chan, V. A. Aksyuk, R. N. Kliman, D. J. Bishop and F. Capasso, Science 291,(2001) 1942. 

[6] U. Mohideen and Anushree Roy, Phys. Rev. Lett. 81, (1998) 4549. 

[7] B. W. Harris, F. Chen, and U. Mohideen, Phys. Rev. A 62, (2000) 052109. 

[8] H. B. Chan, V. A. Aksyuk, R. N. Kliman, D. J. Bishop and F. Capasso, Science 291, (2001) 1942. 

[9] F. M. Serry, D. Walliser and G. J. Maclay, J. Appl. Phys. 84, (1998) 2501. 
[10] R. Esquivel-Sirvent, C. Villarreal y G. H. Cocoletzi, Phys. Rev. A, 64, (2001) 052108 . 
[11] C. Villarreal, R. Esquivel-Sirvent and G. H. Cocoletzi, Int. J. Mod. Phys. A (in press). 
[12] R. Esquivel-Sirvent, C. Villarreal, G. H Cocoletzi and W. L. Mochah, Phys. Stat. Sol. (in press). 

[13] M. Bordag, B.Geyer, G.L. Klimchitskaya, and V.M. Mostepanenko, Phys. Rev. Lett. 85, (2000) 503; GL. Klimchitskaya, 

and V.M. Mostepanenko, Phys. Rev. 63, (2001) 062108. 
[14] E. M. Lifshitz, Sov. Phys. JETP 2, (1956) 73. 

[15] Yu. S. Barash and V.L. Ginzburg, Sov. Phys.-Usp., 18,(1975) 305. 

[16] N. G. Van Kampen, B. R. A. Nijboer, and K. Schram, Phys. Lett. 26A, (1968) 307. 

[17] P. Candelas, Ann. Phys. (N.Y.) 143, (1982) 241. 

[18] R. Matloob, Phys. Rev. A, 60, (1999) 50. 

[19] D. Kupiszewska and J. Mostowski, Phys. Rev. A, 41, (1990) 4636. 
[20] D. Kupiszewska, Phys. Rev. A, 46, (1992) 2286. 

[21] R. Matloob, A. Keshavaraz, and D. Sedighi, Phys. Rev. A, 60, (1999) 3410. 
[22] E. I. Kats, Sov. Phys. JETP 46, (1977) 109. 

[23] V. M. Mostepanenko and N. N. Trunov, Sov. J. Nucl. Physics, 42, (1985) 818. 

[24] V. B. Bezerra, G. L. Klimchitskaya and C. Romero, Phys. Rev. A, 65, (2001) 012111. 

[25] See the review books Photonic Probes of Surfaces, edited by P. Halevi (Elsevier, Amsterdam, 1995) and Spatial Dispersion 
in Solids and Plasmas, Electromagnetic Waves Vol. 1, ed. by P. Halevi (North-Holland, Amsterdam, 1992). For an example 
of a nontrivial application of the surface impedance to non-local excitonic semiconductor superlattices see G. H. Cocoletzi 
and W. Luis Mochan, Phys. Rev. B 39, 8403 (1989). 

[26] M. T. Jaekel and S. Reynaud, J. Phys. I France, 1, 1395 (1991). R. Matloob, A. Keshavaraz, and D. Sedighi, Phys. Rev. 
A, 60, (1999) 3410 . 

[27] G. Plunien, B. Midler and W. Greiner, Phys. Rep. 134, (1986) 87. 

[28] A. E. Gonzalez, Physica A 131, (1985) 228. 



4 



